Renewable and Sustainable Energy Reviews 16 (2012) 2207-2216 


Goes 


ELSEVIER 


Contents lists available at SciVerse ScienceDirect 


Renewable and Sustainable Energy Reviews 


journal homepage: www.elsevier.com/locate/rser 


Towards a de-carbonized energy system in north-eastern Morocco: 


Prospective geothermal resource 


Abdelkrim Rimi &*, Yassine Zarhloule?, Alae Eddine Barkaoui”, Antonio Correia‘, 
Julio Carneiro‘, Massimo Verdoya®, Francis Lucazeau® 


a Institut Scientifique, Université Mohammed V-Agdal, Avenue Ibn Battouta, B.P. 703 Rabat, 10106 Morocco 


b Université Mohammed 1er, Oujda, Morocco 
© Geophysical Centre of Evora (CGE), University of Evora, Portugal 


4 Dipartimento per lo Studio del Territorio e delle sue Risorse, Universita di Genova, Italy 


© Institut de Physique du Globe de Paris, France 


ARTICLE INFO ABSTRACT 


Article history: 

Received 22 April 2011 

Accepted 19 December 2011 
Available online 18 February 2012 


Keywords: 
Geothermal gradient 
Heat flow 

Thermal spring 
Geothermometer 
Geothermal potential 
Morocco 


Geothermal data has been indicating promising potentialities in the north-eastern Morocco. This paper 
presents new temperature data, recently recorded in water boreholes located in the Berkane and Oujda 
areas. Generally, the observed temperature gradients are rather high. One hole near Berkane, revealed an 
average geothermal gradient of more than 110°C/km at depths greater than 300 m. This result confirms 
the geothermal gradient estimated in a mining borehole located about 30 km west of the Berkane bore- 
hole, in which water temperature of 96°C is reached at a depth of about 700 m. Such a high geothermal 
gradient, exceeding by far the ones already determined for northeastern Morocco, could act as a stimulus 
to programs aimed at the geothermal exploitation of high temperature aquifers. 


© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Since 2009, Morocco has been engaged in a sustainable develop- 
ment and renewable energy strategy. With no oil reserves Morocco 
has the necessary assets to achieve these goals, particularly of 
developing free carbon energies, mainly solar, wind and geother- 
mal. Owing to its special geographical and geological position, 
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North East part of Morocco, which is endowed by a natural bounty 
of sunshine and geothermal resources, has seen a number of start- 
ing projects including that of the solar thermal plant of Ain Bni 
Mathar. This paper we will focus particularly on its geothermal 
potential, which may reinforce its development as an emerging 
economic region. 

Since Rimi and Lucazeau [1] have been pointed out the regional 
thermal phenomena (volcanism, hot springs, extensive tectonics 
and seismic activity) in North-eastern Morocco, most advanced 
temperature and geochemical sampling have been done in this part 
of Morocco [2,3]. Work accomplished in this region has highlighted 
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the opportunity for large-scale geothermal anomaly from deep [4], 
while several southern Morocco sedimentary basins such as the 
Souss Coastal Basins are now being explored for lower enthalpy. 
Both subsectors have been inferred to contain enormous quanti- 
ties of stored heat [5,6]. Due to the large number and size of surface 
thermal manifestations (hot springs, recent volcanism and neotec- 
tonic activity), north-eastern Morocco is a region that has attracted 
the attention of geothermal experts since the 1960s [7-9]. 


2. Geological background 


Eastern Morocco is a system of plains, plateaus and mountain 
ranges extending on the eastern parts of the Meseta-Atlas and 
Rif domains. It is represented by the Eastern Middle Atlas that 
continues through the Guercif basin to Beni Snassen, the valley 
of Moulouya, the Oujda-Taza corridor, the “Pays des Horsts” and 
a crustal block without notable alpine deformation, which con- 
stitutes the High Plateaus. Three ENE-WSW structural units are 
defined in Eastern Morocco region Fig. 1: 


- The eastern-rifian foreland composed of Gareb chaotic unit near 
the active Rif chain and a little deformed tectonic foreland (Terni 
Masgout, Beni Bou Yahi, Beni Mamou and Beni Snassen), pro- 
longed towards the east (in Algeria) by the Traras mounts. 

- The Guercif and Taourirt-Oujda basins which show sediments of 
middle Miocene to Quaternary age, and are connected to the Tafna 
basin to the east (Algeria). 

- The Taourirt-Oujda mounts which show sub-horizontal Juras- 
sic sediments affected by fault systems whose development and 
geometry are strongly dependent on the pre-existing structures 
in hercynian basement. They are prolonged to the east by the 
pre-tellian foreland of Tlemcen. 
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The Rif domain, located northwards, belongs to the Betic-Rif belt 
and is part of the peri-Mediterranean alpine chain. With the Atlas, 
these morpho-structural zones result from the superposition of 
several tectonic, compressive and extensional phases of alpine age. 

These phases have occurred since the Late Cretaceous. Most of 
these zones are ENE-WSW to E-W directed (e.g., “Pays des Horsts”) 
and NE-SW (Middle Atlas, valley of Moulouya), reflecting part of the 
alpine geological evolution. 

From Landsat image analysis, El Hammichi et al. [10] had 
depicted four faults systems: 


(a) The first system, NE-SW oriented, is well developed and rep- 
resented by several faults with regional importance. Along 
this fault, structures are very clear and indicate prevailing left 
strike-slip movement; 

(b) NW-SE (N135°E to N160°E) trending faults are often longer 
than 15 km and extend the southern sector; 

(c) The third system in this sector is represented by arcuate 
ENE-WSW to E-W trending faults; 

(d) The last fault system, affecting especially the western side of 
the studied area, is represented by faults trending N-S. These 
faults are clear and long lineaments. Kinematic indicators and 
vertical offsets testify normal movements. 


2.1. Recent volcanism 


On the scale of Morocco, the volcanic centers extend roughly 
along a NE-SW direction since the Mediterranean coast to the Anti- 
Atlas and the Canary Islands. 

In North eastern Morocco, late Miocene to Pliocene volcanic 
activities correspond to a shoshonitic edifices in the Guilliz and 


200000 


N 
À ALBORAN SEA 


Le 


£31) 


3 


SX 


asto fit // 


D” | 


M Foriand (African plate) e  60<HFD<100 


EM Ritan Thrust units @ :00<HFD<140 
Post nappes(neogene and Quaternary) 

HEEB Miocene and Plio-Quatemary 8 140<HFD 

Fault 


Fig. 1. Simplified structural sketch map and volcanism in northeastern Morocco (modified after Hernandez [29]; Hervouet [30]). Heat flow density values are superimposed 
on this map; numbers in circle represent spring temperature values, K: mining borehole at Kariat Arekman, B: borehole at Berkane and F: borehole at Fezouane site. 
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Koudiat Hamra (near Taourirt) (Fig. 1) with age ranging from 8 to 
4.5 My [11] and calcalkaline centers in the Oujda area (Fig. 1) with 
age ranging from 6.2 to 4.8 My [12]. 

Regionally, the study area is located east of Gibraltar Arc, which 
is currently the subject of many geological, geophysical and geo- 
dynamical investigations. This has shown that the Gibraltar Arc is 
characterized by thinned crust and lithosphere, with a fast thick- 
ening towards the Betics and the Rif mountain ranges. 

In eastern Morocco, magmatism developed during the Miocene. 
A more recent episode consisting of alkaline volcanic activity has 
taken place from 6 to 0.8 Ma, and occurred along a NE-trending 
belt extending from the Trans-Alboran zone to the Middle Atlas 
domain and to the Anti-Atlas (J. Siroua volcano). Their geochemi- 
cal signature is similar to that of intra-oceanic island basalts, and 
suggests the role of an asthenospheric “hot lineament” [13]. Such 
interpretation is consistent with gravimetric and geodetic model- 
ing of the lithosphere [14-16]. This MHL could extend from the 
Canary Islands to southeast Spain at least. 

Lithospheric thinning under the Alboran Sea and thickening 
under two ridges of the adjacent Gibraltar Arc is of basic impor- 
tance for geodynamic models which invoke convective removal 
and orogenetic collapse, a mantellic delamination and an active 
subduction [17-21]. The Alboran Sea continues eastward with the 
Algerian oceanic crust basin [22]. 


3. Hydrogeological and geochemical features 
A sub-surface reservoir is suitable for extracting heat, if its 
porosity, permeability, thickness, lateral extension and tempera- 


tures are sufficiently large. For a preliminary description of the 
geothermal areas of Morocco, two criteria are considered: 
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- The structural, lithological and hydrogeological context to estab- 
lish the existence or the absence of aquifers; 

- Temperatures of aquifers, or impermeable formations up to 3 km 
deep. 


The variety of tectonic and lithostratigraphic in Morocco has 
created various hydrogeothermal systems with heterogeneous 
lithological and physical properties. Regional studies in Prerif [23], 
eastern Morocco [24], the plateau of phosphates [25], Essaouira- 
Agadir Basin [26,27] and the central Morocco [28] have been 
gathering background information on aquifers from archives of 
piezometers and oil Boreholes: structure, depth, porosity, modflow 
and groundwater velocity. 

But features such as permeability, transmissivity and salinity are 
often lacking. According to these studies, the potential geothermal 
aquifers resources in Morocco are of two types: 


(1) In the latest Tertiary and Quaternary deposits, water pumped 
out from pore spaces during compaction could constitute up to 
700 m depth, aquifers of good low temperature (50-60 °C). 

(2) The typical system of water circulation for pre-Tertiary aquifers, 
1000-4000 m depth, is guided by differences in hydraulic load 
associated with formation changes. The meteoric water infil- 
trates into the hills and topographic relief and discharges in 
areas of low groundwater level, mainly in Mesozoic carbon- 
ates of the major basins. If the permeabilities are favorable, the 
water can penetrate deep enough and warm up gradually (aver- 
age temperature is 100°C). However, the lithological nature 
may lead to the formation of local or even regional convective 
systems. The descent of cold water and the rise of hot water can 


disturb the thermal field. 
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Fig. 2. Oujda and Berkane Liassic aquifer [4]. 
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Table 1 

Thermal data from hydrogeological and petroleum boreholes. 
Descriptives codes Latitude Longitude Elevation Maximum depth Geothermal gradient Heat flow 

(m) (m) (°C/Km) (mW/m?2) 

Fezouane2 34.927 —2.213 252 49.0 140 
1624/7 Berkane 34.93 —2.3 160 127.0 510° 
2952/12 34.693 —1.818 560 500 49.0 111 
2362/12 34.647 —1.930 215 43.0 87 
1199/12 34.645 —1.964 641 125 34.0 90 
1225/12 34.782 —1.930 75 49.0 98 
2899/12 34.642 —1.847 620 128 14.0 89 
186/16 34.015 —1.913 385 255 32.0 60 
A4 34.305 —2.227 1068 260 39.0 81 
A13 34.311 —2.232 1095 170 30.0 119 
A 14 34.309 —2.238 1068 350 34.0 100 
A7 34.299 —2.222 216 36.0 89 
P3 34.308 —2.232 240 33.0 68 
BR41W 34.45 —1.802 1058 170 45.0 144 
TEFADIS 34.507 —1.780 1125 100 45.0 112 
BN38NW 34.452 —1.81 1053 100 43.0 105 
BN30 34.443 —1.792 1050 120 37.0 108 
BS38N 34.450 —1.804 120 36.0 84 
BR41S 34.448 —1.807 120 35.0 91 
NADOR1 35.326 —2.662 2703 41.8 86 
AREKMAN 35.109 —2.739 678.7 124.0 192 
YADIM 35.093 —2.475 70 40 57 
Grf1 34.279 —3.424 36.3 92 
Tafrata1x 34.22 —3.106 34.6 76 
Msd1 34.245 —3.47 37.0 90 


* This value is still a rough estimate since we have no real determination of thermal conductivity and porosity. Ranging from 3.7 to 5.5 W/m°C, the thermal conductivity 


of the dolomite decreases as the porosity and temperature increase. 


For an area to be prospective for high enthalpy geothermal 
resource there must exist high average geothermal gradients to 
create high temperatures at depths over a significant area which 
can be economically drilled. In addition to this the heat must be 
able to be economically exploited. This requires the existence of 
suitable reservoir rocks with either high natural porosity and per- 
meability or the conditions that allow permeability to be artificially 
enhanced by reservoir engineering. 

The sedimentary formations in north-eastern Morocco units 
contain hot reservoirs; in particular Liassic carbonates up to 500m 
thick which is the most important aquifer in the region (Fig. 2). The 
main aquifer is fed into the Horsts area and discharges in the plain of 
Oujda-Angad, and is composed of carbonate rocks (limestones and 
dolomites) from the lower Cretaceous, the Lias (Fig. 2). The isobaths 
of Lias vary from a few meters south to 1000 m towards the north. 
Hydrostratigraphic studies based on geological and hydrogeologi- 
cal data from boreholes show the Liassic carbonate to be the main 
hydrogeothermal reservoir of the hot springs in the region. 

The Liassic reservoir of Eastern Morocco belongs to the Atlas 
domain, which is characterized by a succession of NE-SW horst 
and graben of pre-Miocene age and by Plio-quaternary basaltic 
volcanics, associated with crustal thinning (Table 1). 

This reservoir is at the origin of twenty-four thermal springs. 
Table 2 whose temperature ranges from 20 to 54°C and the flow 
rate can reach 401/s at certain points. Salinity varies from 0.28 to 
29.6 g/l. Some of these hot springs play an important role in the 
economy of the area (e.g., Fezouane, near Berkane and Hammam 
Ben Kachour, at Oujda [3,4]. 

Fig. 2 shows the Liassic aquifer in both Oujda and Berkane region. 
The piezometrical map shows that flow direction is from SE to NW. 
The recharge area of the confined aquifer is localised in the south, 
where the Liassic outcrops. In the Fezouane region, the aquifer is 
artesian. 

In the Oujda region, the Liassic limestone forms the main aquifer. 
Limestone outcrops in Beni Snassen, High-Atlas, Middle-Atlas and 
the southern part of the basin. The thickness varies from 50m to 
1140 m, and the maximum depth is 1370 m. The groundwater flow 
is from south to north of basin Fig. 2. The aquifer temperature and 


salinity vary from 20°C to 52°C and from 130 mg/l to 3000 mg/l, 
respectively. 

According to Zarhloule [3], the hot temperature and the arte- 
sian rise of most of the thermal springs in north-eastern Morocco 
originate from deep circulation in relation to a system of basement 
active faults, individuating horsts and grabens, in a recent volcanic 
area. 

Winckel [28] performed a thorough geochemical analysis of the 
main hot springs in Morocco and found that eleven of those located 
along the MHL (Moroccan Hot Line) release CO: that is partially 
of deep origin. Recently, Tassi et al. [31] confirmed that CO2-rich 
thermal springs with 3He anomalies are likely related to this hot 
line. They are mainly distributed along a NE-SW trend from Nador 
to Taza, and from Fes (Moulay Yacoub) to Oulmes south of the 
Rif frontal thrust. The contemporary presence of 3He anomalies 
and minor recent basalt outcrops indicate that CO; originates from 
mantle degassing or deep hydrothermal systems in these thermal 
discharges 

The geochemical analyses and water classification by means of 
the diagram of Langelier and Ludwig gave the following hydrogeo- 
chemical facies Ca (Mg)—HCO3), NaCl and Ca-Mg [6]. 

The estimated temperatures by silica geothermometers vary 
between 37 and 85°C (quartz of Fournier), 23 and 73°C (quartz 
of Arnorsson), 54°C (chalcedony of Fournier) and 57°C (chal- 
cedony of Arnorsson). The measured emergence temperatures vary 
between 20 and 54°C. The correction of the water estimated reser- 
voir temperatures by the mixture geothermometer reveals that the 
temperature in the original reservoir may exceed 110°C, witha rel- 
atively low rate of mixture between hot deep water and cold surface 
water. 


4. Geothermal pattern and new borehole thermal logs 


Hydrogeothermal structure was impulsed by a geothermal 
anomaly measured in a mining borehole drilled in rhyolitic and 
andesitic series at Kariat Arekman (35.11°N, 2.74°W; Fig. 3). The 
well, 680 m deep, revealed from 450m depth salt artesian water 
(20-25 g/l and a flow rate of 1.4-21/s) with a temperature of 
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Table 2 
Thermal springs properties. 
Name Latitude Longitude  T(°c) Rate Gas Chimical Facies Reservoir Use Project 
(I/s) 
Sidi chaffi 34.9 —3.61 32 4 CO2 Na-HCO3 Marl-limestone Drink, therapeutic Hammam 
Midar cretaceous bathing, irrigation 
Ain Messouda 35.03 —3.25 17-21 3 Chloride, sodium, Gastric and urinary  Botling 
sulfate, diseases manufactory 
calcic-magnesian project in 1959 
Ain Hammam 35.17 —3.19 32-35 1 H2S, CO2  Na-C Miocene rhyolite Dermatologica Hammam 
Beni Sidel post nappe, Tuffset bathings 
breccias 
interbedded 
Hammam 35.17 —3.09 34 2 Na-C Massive limestone Dermatologica Hammam 
Chaabi and gray slate of bathings 
Jurassic of 
Taliouîne 
Hajra Safra 34.26 —3.56 2 CO2 Na-C Marl-limestone Dermatologica Hammam 
cretaceous bathings 
miocene 
Ain Haddou Ou 35.14 —3.09 29 1 H2S, CO2 Na-C Massive limestone Dermatologica Hammam 
Hamma and gray slate of bathings 
Jurassic of 
Taliouîne 
Ain Goutitir 34.35 —3.06 44-49 12 CI-Ca Limestone- Hydrotherapy Hammam 
dolomitic lower 
lias-Domerian 
Gufait Springs 34.24 —2.38 26 10-15 SO4-Ca Limestone- Bathing, drink and 
dolomitic irrigation 
Lias 
Ain Bourachad 33.86 —3.60 32 2 Na-Cl Limestone- Dermatology Two pools 
dolomitic 
Lias-Dogger 
Ain El Hamra 34.61 —3.91 24-26 6-15 CO2 HCO3-Ca Marl-limestone Repute for kidney Hammam 


and sandstone diseases 


Miocene 


42°C at the surface, and more than 90 °C at the bottom [32]. This 
good result infirmed the earlier studies [8,9], which concluded 
on the base on 60 meters deep thermal profiles, that north- 
eastern Morocco does not constitute really a zone with significant 
geothermal potential. Carried out in the very shallow boreholes, 
such works cannot be considered conclusive, because the ther- 
mal state is strongly disturbed by water movement at this depth. 
In addition, more recent deep geophysical investigations (seismo- 
logical, gravity and magnetic) claim that the positive geothermal 
anomaly in the Rif and Eastern Morocco is explainable as a surface 
manifestation of a deep and wide tectono-magmatic lithospheric 
process [33]. 

The surface heat flow and broad thermal conductivity structure 
(and thus the temperature field) of northeastern Morocco has been 
estimated since 1987, [1,3,4,6] provided broad approximations of 
northeastern Morocco continent temperature field, but suffered 
from a sparse geographic distribution of heat flow and thermal 
conductivity data. 

The analysis of the long wavelength anomalies highlights a 
heat flux increasing northeastward: 80-110 mW/m? in the east- 
ern Rif, northeastern Morocco, Alboran Sea (east of the meridian 
line 4°W), southeastern Spain and northwestern Algeria. In east- 
ern Morocco, the crust is thinned to 15-20km. The highest heat 
flow (80-140 mW/m2) and geothermal gradient (35-50 °C/km) are 
observed in the north-eastern part of the country (Middle Atlas and 
Eastern Rif), which is characterized by intensive Quaternary volcan- 
ism, high residual magnetic anomalies and shallow Curie depths 
[34]. The Gibraltar Arc region is characterized by a radial heat flow 
pattern with increasing values from the outer ranges towards the 
center and eastern part of the basin [2,35] (Fig. 4). The most strik- 
ing geothermal feature is the abnormally hot zone observed east of 
4°W and the heat flow increase towards the Algerian basin, where 
we find volcanic environments and high values typical of a young 
oceanic lithosphere [2,36]. 


Fig. 5 shows a set of temperature data recorded in boreholes of 
the study area in 1986 [1]. It shows that the geothermal gradient is 
generally larger than 30°C/km. New thermal logs were performed 
in 2007 and 2009 at some boreholes (see Fig. 1 for their location). 
Fig. 6 points out abnormally hot deep geotherms beneath Eastern 
Rif and Alboran basin [37]. 

With the objective of improving the geothermal dataset and 
gaining a better insight into the geothermal potential of eastern 
Morocco, an exploratory geothermal survey was conducted in 2007 
and 2009; in the region of Oujda (Fig. 6) to log some hydraulic 
boreholes. In particular, two wells in the region of Berkane show 
high geothermal gradients and could prove valid for high enthalpy 
geothermal purposes. This paper is a summary of the results of 
these recent investigations and discusses the new data in the light 
of previously collected information in deep and shallow boreholes 
(Fig. 7). 

The borehole 2952/12. Fig. 8a shows a good concordance 
between the thermal profiles measured after 20 years. At 600 m 
depth (corresponding to the lithological change from marl to car- 
bonate), the borehole may be invaded with by a warm water arrival. 
This results in a sharp decreases and even an inversion of the 
geothermal gradient. The thermal log of the well located 30 km west 
of Berkane. Fig. 8b is characterized by an increase of the geother- 
mal gradient at 300 m depth. Notice that the increase of geothermal 
gradient cannot, to a certain extent, be explained by a change in 
lithology. In that case, as dolomite should have thermal conduc- 
tivity higher than clay, the geothermal gradient should decrease. 
For the moment, no thermal conductivity measurement is available 
and no heat flow density estimate can be done. 

The thermal log of the Fezouane hole (Fig. 8c) exhibits a ther- 
mal gradient, which apparently varied with depth. The thermal 
gradient is 42 °C/km in the uppermost 60-120 m section, and then 
decreases to 14°C/km between 120 and 150 m, to increase again 
to about 70 °C/km in the deeper part of the hole. Notice that this 
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Fig. 3. Observed temperature in the Kariat Arekman hole [32]. Local estimated heat flow density exceeding 200 mW/m? [1]. 


log is significantly affected by water flowing within the hole, since 
it crosses a confined, artesian aquifer. 


Heat Flow Density Variations in Morocco 
and neighbouring areas 


4.1. Heat flow 
35 If we assume that the conduction is the main mode of heat trans- 
fer in the crust, average vertical heat flow, at the earth’s surface can 
be determined by: 


Qs = Azgrad, + 5 H 
30 
where (z is the thermal conductivity (W/mK) at depth z and gradz 
the geothermal gradient (K/m) at depth z, (His the additive radioac- 
tive heat production between the surface and the depth z. 


4.2. Some results for the Berkane well 


This change in the geothermal gradient is not a result of a 
thermal conductivity change. A simple calculation shows that 
the average geothermal gradient above 300 m is about 29°C/km 
Mauritania (0.0293 °C/m) while the geothermal gradient below 300 m is about 
127°C/km (0.1275°C/m). At a depth of about 470 m the tempera- 
ture is 50°C, and assuming that the same formation reaches depths 
Fig. 4. Terrestrial heat flow density in Morocco and neighboring regions [34]. of 700 m (see lithologic log) this means that an estimated temper- 
ature of about 78°C is attained at about 700 m. 

If we assume a thermal conductivity of 1.5 W/mK for the clay 
above 300m and a thermal conductivity of 4.0W/mK for the 
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Fig. 5. Thermal profiles in mining and water boreholes done in Oujda-Berkane 
region. 


dolomites below 300 m we can conclude that the heat flow density 
(HFD) along the well is not constant and therefore there is geother- 
mal energy trapped for levels below 300m deep. The estimated 
HED above 300 m is 43 mW/m2, and below 300 m is 510 mW/m2. 

If we assume that the well is representative of the geologic, 
hydrogeologic, and thermal status for a relatively large are (let’s 
say tens of square km), we can try to estimate the heat in place (Ho) 
within the aquifer using the following equation: 


Ho = [C1 — P)0mCm + PewCw] - (Tr — To) -A - Az 


Temperature (°C) 
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Fig. 6. Crustal geotherms beneath Eastern Rif (RF) and AlBoran basin (AB) [37]. 
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Fig. 7. Geothermal gradient (a) and Heat flow density (b) distribution in the study 
area. 


@ is the effective porosity, C is the mean specific capacity (j/KgK), 
p is the mean density of the rock column (kg/m3), T; is the tem- 
perature at the top of the aquifer (°C), Tọ the average regional 
temperature at the surface of the earth (°C), a is the area of the 
aquifer, and À, is the thickness of the aquifer; the subscripts m 
and w refer to rock matrix and water, respectively. 

For the Berkane well if we assume that: ®=0.1, Tp =20°C, 
Tr=30°C, pw=1000kg/m?, pm=2600kg/m?, Cy =4190J/kgK, 
Cm = 840J/kgkm the heat in place per square meter is about 
9.5 x 109 J/m? (9.5 GJ/m?). 

It must be said that we were able to go with the temperature 
datalogger as deep as 900m in the well; however, we were not 
counting that such high temperatures were reached at depths of 
about 500m and so the temperature data logger saturated just 
above 50°C. The next temperature log must then be obtained with 
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Fig. 8. (a) Thermal profile measured in the well 2952/12 (34.4 N, 1.5°W) in 1986 and 2007. (b) Temperature log of the well 1624-7 in the region of Berkane. (c) Temperature 


log in the well Fezouane2 in the region of Berkane. 


a temperature data logger that allows measuring temperatures 
higher than 50°C or with another temperature measuring device 
that can measure temperatures higher than 50°C. 


4.3. Present day thermal waters direct use in northeastern 
Morocco 


Tough geothermal energy research in Moroccois relatively more 
than forty years old [38], the utilization of geothermal resources 
has not yet been enough developed and still restricted to sparse 
low enthalpy direct uses. Among the geothermal potentialities of 
Morocco, numerous hot springs occur and are mostly concentrated 
in the northern country: Central Morocco, Prerif, and northeastern 
Morocco where high geothermal gradients are observed. Because 
of the low to intermediate enthalpy level of these geothermal hot 
springs, the direct use of the warm springs is still for therapeutic 
or for their chemical contents. However except for the spring of 


Moulay Yacoub of Fez where modern therapeutic bathing exists, the 
thermal waters in Morocco remain still primitively used. In spite of 
the whole energy dependence of Morocco and the development of 
the geothermal studies, the Moroccan authorities are not decided 
yet to promote the use of the geothermal energy contained in the 
underground of the country which could make primarily profit it’s 
agricultural and fishing vocation. 

For a flow rate of 51/s and an average of 31/s and a temperature 
change of 10° C and a total of 12 natural spring baths (Benkachour, 
My Yacoub, etc.). 


Installed capacity (MWt) = peak use x “AT” x 0.004184 
x number of baths = 101/s x 10°C x 0.004184 x 12 = 2.51 MWt 


Annual energy use (TJ/year) = average user“ AT” x 0.1319 
x number of baths = 31/s x 10°C x 0.1319 x 12 = 47.48 TJ/year 
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Therefore geothermal annual used only by the natural hot baths 
in Morocco, remains less than 50TJ/year which corresponds to a 
capacity of 3MWt. 


5. Discussion 


The new thermal logs recently recorded show possible ther- 
mal perturbation due to groundwater effects and/or water flowing 
within the holes. The acquired data need further processing and 
careful analyses before their use for geothermal modeling. How- 
ever, some general conclusions can be drawn from the observed 
temperatures. 

A change in the geothermal gradient does not appear as a result 
of a change in thermal conductivity. A simple calculation for the 
Berkane well. Fig. 8b shows that the average geothermal gradi- 
ent above 300 m is about 29°C/km while the geothermal gradient 
below 300 mis about 127 °C/km. At a depth of about 470 m the tem- 
perature is 50°C, and assuming that the same formation reaches 
depths of 700 m (see lithologic log) this means that an estimated 
temperature of about 78°C is attained at about 700 m. 

If we assume a thermal conductivity of 1.5 W/mK for the clay 
above 300m and a thermal conductivity of 40W/mK for the 
dolomites below 300 m we can conclude that the heat flow density 
along the well is not constant and therefore there is heat trapped 
at depth below 300 m. 

The hot water may be recovered essentially in Lias series, 
for geothermal purposes. However, the problem with the Liasic 
aquifers is that they are formed by intrinsically compact limestones 
and dolomites, but which can contain water in fractures and even of 
karsts systems. The formations may be thus transmissives but not 
very accumulatives which makes hardly the definition of the thick- 
ness of the aquifer, the porosity and the transmissibility. Further 
investigations must be deserved for a good estimation of the thick- 
ness and the porosity of the aquifer which remain key parameters 
for any attempt of geothermal spring assessment [39,40]. 

The borehole of Berkane may not be considered as represen- 
tative of the geologic, hydrogeologic, and thermal status for a 
relatively large area, owing probably to complex upward or down- 
ward groundwater flow. We remark the difference in thermal 
gradient between the two boreholes, Berkane (Fig. 6) and Fezouane 
(Fig. 8c), between 150 and 200 m depth. 

Notice that we were able to go with the temperature data logger 
as deep as 900 m in the well; however, at depths of about 500 m and 
so the temperature data logger has reached its maximum operation 
limit (50°C). Future investigations should then use a temperature 
data logger that allows to measure temperatures higher than 50°C. 


6. Conclusion 


It was hypothesized from geological and geophysical thermal 
spring data that Northeastern Morocco is likely to be a region 
of high average geothermal gradients and thus prospective for 
geothermal energy. 

The heat flow measurements in hydraulic, mining and oil bore- 
holes returned results that have validated this hypothesis. 

The sediments in the study area have moderate thermal conduc- 
tivity, and have been sampled for some laboratory measurements 
and estimated by geometric models to be between 2.5 and 7 km in 
total thickness. With this amount of moderate conductivity sed- 
iments overlying a thinned crust, plus its proximity to existing 
power infrastructure, the Northeastern Morocco is now considered 
prospective for economic geothermal energy development. 

The study area is part of a larger, elongate, geographical region 
extending from Morocco to Alboran Sea and South-eastern Spain. 

These large resources of geothermal energy in northeastern 
Morocco would not be unnoticed. It can be believed they will attract 


interest of potential investors and will be profitable also for those 
whole areas. 

High geothermal gradients are observed in NE Morocco In the 
region of Berkane, temperatures higher than 50°C were measured 
in a water well at depths of about 470m while this temperature 
was almost reached at 600 m in another well of Oujda region. 

This large geothermal anomalies occurring at depth and this 
already observed in a mining borehole at Arekman area (Eastern 
Rif) indicate that medium to high enthalpy uses are not discarded. 

Other wells are available in northeastern Morocco for tempera- 
ture logging, so that the geothermal potential can be estimated. 
For a valuable assessment of the potential geothermal springs, 
detailed geophysical surveys have to be carried out in the 
Oujda-Berkane-Arekman area. 
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